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Non-Singlet QCD Analysis of the Structure Function in 3-Loops 
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First results of a non-singlet QCD analysis of the structure function Fi{x, Q 2 ) in 3-loop order based on the non- 
singlet world data are presented. Correlated errors are determined and their propagation through the evolution 
equations is performed analytically. The value for a s (Mz) is determined to be 0.1135 + / — 0.0023/0.0026 
compatible with results from other QCD analyses. Low moments for u v (x), d v (x) and u v (x) — d v (x) with correlated 
errors are calculated which may be compared with results from lattice simulations. 



1. INTRODUCTION 

A consistent 3-loop QCD analysis of the unpolar- 
ized structure function F 2 (x,Q 2 ) can be carried 
out after having the recently completed next-to- 
next-to-leading order (NNLO) anomalous dimen- 
sions available 1 in addition to the 2-loop Wil- 
son coefficients Two NNLO QCD analyses 
have been performed previously based on struc- 
ture function and other hard scattering data pj], 
and structure function data only @|. In these 
analyses singlet- and non-singlet evolution was 
dealt with in parallel. The present analysis con- 
centrates on the non-singlet evolution only in or- 
der to firstly obtain an accurate as possible pic- 
ture for the valence quark distributions and the 
value of the QCD coupling constant a s (Af|) with 
correlated errors in NNLO. Analyzing only the 
non-singlet data has the advantage that large 
gluon and sea quark effects remain decoupled. A 
NNLO QCD analysis of F 2 (x,Q 2 ) allows to re- 
duce the theoretical error in determining a s (Q 2 ) 
to at least the level of the experimental error since 
the factorization and renormalization scale uncer- 
tainties reduce significantly. Comparison of QCD 
analysis results with results from recent lattice 
simulations concerning the low order moments 
has shown astonishing agreement in the polar- 
ized case 0. With the steadily improving lat- 
tice calculations this might also become feasible 



in the unpolarized case, where systemantic effects 
have still been large during the last years. In 
this letter we describe the results of an analysis 
of the deeply inelastic non-singlet world data for 
charged lepton-nucleon scattering. 

2. QCD FORMALISM 

In Mellin-AT space the non-singlet (NS) parts of 
a structure function Fi(N, Q 2 ) are given by 



F^{N,Q 2 ) 



[l + Cu(N)a + C 2l (N)a 2 ] 
/ ±lV (A,Q 2 ), (1) 



where the Cki(N) are the corresponding Wil- 
son coefficients and f ±,v (N,Q 2 ) are the non- 
singlet quark combinations. The symbol a de- 
notes the strong coupling constant normalized to 
a(Q 2 ) = a s (Q 2 )/4:ir. The quark combinations to 
be considered are 

A ± = (u±u) - (d±d) , (2) 
v = (u - u) + (d - d) . (3) 
The non-singlet parts of F 2 are proportional to 



F 2 NS oc V, 
2 18 



and 



(4) 
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F^ s stands for the difference of proton and 
deuteron data in the range x < 0.3 while the 
other combinations are used in the valence ap- 
proximation for x > 0.3. All these combina- 
tions evolve as -(—combinations in Q 2 . The rele- 
vant parton densities to be determined in a non- 
singlet QCD analysis are xu v (x,Q 2 ), xd v (x,Q 2 ), 
and x{d — u){x, Q 2 )- 

The solution of the evolution equation to 3- 
loops reads for the valence distribution as a ex- 
ample 
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where the Pi are the splitting functions in Mcllin 
space. 

3. PARAMETERIZATION 




Figure 1. The parton density xu v at the input 
scale Q 2 — 4.0 GeV 2 (solid line) compared to re- 
sults obtained by MRST (dashed-dotted line) 
and Alekhin (dashed line) @|. The shaded areas 
represent the fully correlated la statistical error 
bands. 



since their respective errors are still large. Fur- 
thermore, the analysis requires different u and d 
quark densities. We adopted the choice 6 which 
gives a good description of the [d — u) data from 
E866 



The parameterizations of the above mentioned 
parton densities at the input scale of Qq = 
4.0 GeV 2 are determined as 

xu v (x, Qq) = % u (1 a; ) u 

(1 - l.lOSa;* + 26.283a;) , (6) 

and 

xd v (x,Q 2 a ) = A dv x ad (l-x) b * 

(l + 0.895a;i + 18.179a;) . (7) 

Here the values for the coefficients in the polyno- 
mial given as numbers in the intermediary range 
of x are obtained by a fit and are then kept fixed 



x{d-u)(x,Q 2 ) = l.igsx 1 - 24 ^-^) 9 - 10 

(1 + 14.05a; - 45. 52x 2 ) . (8) 

The normalization constants A Unj and are 
fixed by the conservation of the number of va- 
lence quarks: J Q u v (x)dx = 2 , J d v (x)dx = 1. 

The remaining four parameters are determined 
in the fit for u v and d Vl the low-a; and high-a; pa- 
rameters a and 6, respectively, along with Aqcd- 

4. DATA 

The results presented here are based on 762 data 
points for the structure function F2(x,Q 2 ) mea- 
sured on proton and deuteron targets. The ex- 
periments contributing to the statistics were : 
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Figure 2. The parton density xd v at the input 
scale Qq with the same conditions as in Figure 
la. 



BCDMS 8 1, SLAC 0, NCM QUI, HI Q3], and 

ZEUS Q21. The BCDMS data were recalculated 
replacing Rqc d by -R1998 . All deuteron data 
were corrected for Fermi motion and offshell ef- 
fects ^31- The non-singlet structure function 
F^ s was constructed according to the relation 
F^ s — 2(i*f - F2) from proton and deuteron 
data given at the same x and Q 2 values. The 
region x > 0.3 which was selected for F,f ' d is ex- 
pected to be dominated by valence quark con- 
tributions while for the region x < 0.3 F^ 8 
data were used. In addition, a Q 2 region of 
4.0 < Q 2 < 30000 GeV 2 was chosen and a cut 
in the hadronic mass W 2 > 12.5 GeV 2 was ap- 
plied in order to diminish higher twist effects. 

In the fitting procedure we allowed for a rel- 
ative normalization shift between the data sets 
within the overall normalization uncertainties 
quoted by the experiments or assumed accord- 
ingly. These normalization shifts were fitted once 
and then kept fixed. 



5. CORRELATED ERROR CALCULA- 
TION 

The systematic errors which are known to be 
partly correlated are not treated separately in this 
analysis. For all data sets we used the simplest 
procedure by adding the statistical and the to- 
tal systematic errors in quadrature being aware 
of the fact that this will eventually overestimate 
the errors. 

The fully correlated ler experimental error 
bands are determined via Gaussian error prop- 
agation demanding that only fits with a positive 
definite covariance matrix are accepted. The gra- 
dients of the parton distributions with respect to 
the variable parameters needed here can be cal- 
culated analytically at the input scale Q\. Their 
values at Q 2 are given by evolution in Mellin-JV 
space. 

6. RESULTS 

The fit results are summarized in table 1. 



u v 


a 
b 


0.314 ± 0.007 
4.199 ± 0.032 


d v 


a 
b 


0.413 ± 0.047 
6.196 ± 0.332 


A (4) 
QCD 


227 ± 30 MeV 


x 2 /ndf = 652/757 = 0.86 



Table 1: Parameter values of the NNLO QCD fit 
based on the non-singlet world data on F 2 em (x, Q 2 ). 

The value of ^q CD is quite stable against a 
variation of the Q 2 cut on the data when varying 
it between 4.0 GeV 2 to 10.0 GeV 2 . From the 
fitted value of ^qcd the following value of a s is 
extracted: 

a s {M%) = 0.1135 («!>*)• 

This value is within the errors in good agreement 
with results from other NLO/NNLO QCD analy- 
ses, see Table 2, and with the latest value for the 
world average of 0.1182 ± 0.0027 \Hj\. 
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The resulting parton densities xu v (x) and xd v (x) 
at the input scale of Qq = 4.0 GeV 2 are pre- 
sented in Figures 1, 2. Comparison with results 
from global analyses shows satisfactory agree- 
ment. While xu v is rather well determined, the 
error band for bit broader as also found 

by the other analyses. 





a s (M'z) cxpt theory 


Rof. 


NLO 






CTEQ6 

MRST03 

A02 

ZEUS 

HI 

BCDMS 


0.1165 ±0.0065 

0.1165 ±0.0020 ±0.0030 

0.1171 ±0.0015 ±0.0033 

0.1166 ±0.0049 

0.1150 ±0.0017 ±0.0050 

0.110 ±0.006 


[17] 
[3] 
[4] 
[19] 
[12] 
[9] 


BB (pol) 


+0 009 

0.113 ±0.004 n„"« 
—0.006 


[5] 


NNLO 






MRST03 
A02 

SYOl(ep) 
SYOl(^N) 


0.1153 ±0.0020 ±0.0030 
0.1143 ±0.0014 ±0.0009 
0.1166 ±0.0013 
0.1153 ±0.0063 


[3] 
[4] 
[19] 
[19] 



Table 2: Comparison of NLO and NNLO results 
on a s (Mz) from deep inelastic scattering, including 
also global fits. The NLO HI value is subject to an 
additional error of +0.0009/ - 0.0005 and the NLO 
ZEUS value of ±0.0018 due to model dependence. 



Low moments for u v , d v and u v — d v with cor- 
related la errors have been calculated. An in- 
teresting comparison can be made for the second 
moment of u v — d v with a result from recent lat- 
tice simulations. A value of 0.180 ± 0.005 was 
derived here and can be compared with results of 
upcoming lattice simulations, cf. e.g. |16| . This 
comparison should be performed for higher mo- 
ments in the future as well as for the moments of 
the individual distributions u v and d v . 

7. CONCLUSIONS 

A non-singlet QCD analysis of the structure 
function F2(x,Q 2 ) based on the non-singlet 
world data has been performed at 3-loop order. 
The value determined for a s (M§) is compatible 
within the errors with results from other QCD 
analyses and with the world average. New pa- 
ramcterizations of the parton densities u v and d v 
including their fully correlated la errors are de- 



rived. Low order moments for u v , d v and u v d v 
are calculated from this analysis. It will be in- 
teresting to compare them with upcoming results 
from lattice simulations. 
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